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ABSTRACT 

The goal of the Herschel Open Time Key programme "TNOs are Cool!" is to derive the physical and thermal properties for a large sample of 
Centaurs and trans-Neptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. We present results for seven 
targets either observed in PACS point-source, or in mini scan-map mode. Spitzer-MlPS observations were included for three objects. The sizes of 
these targets range from 100 km to almost 1000 km, five have low geometric albedos below 10%, (145480) 2005 TB190 has a higher albedo above 
15%. Classical thermal models driven by an intermediate beaming factor of 77=1.2 or 77-values adjusted to the observed colour temperature fit the 
multi-band observations well in most cases. More sophisticated thermophysical models give very similar diameter and albedo values for thermal 
inertias in the range 0-25 J nr 2 s~ 05 K _1 , consistent with very low heat conductivities at temperatures far away from the Sun. The early experience 
with observing and model strategies will allow us to derive physical and thermal properties for our complete Herschel TNO sample of 140 targets 
as a benchmark for understanding the solar system debris disk, and extra-solar ones as well. 
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1. Introduction 

Trans-Neptunian objects (TNOs) are believed to represent one of 
the most primordial populations in the solar system (Morbidelli 
et al. 2008 '). The TNO population comprises (i) the main Kuiper 
belt beyond the orbit of Neptune (~32 - 50 AU), consisting of 
objects in resonant and non-resonant orbits, and (ii) the halo out- 
skirts of "scattered" and "detached" bodies beyond 50 AU. The 
Centaurs, an unstable orbital class of minor planets (e.g., Horner 
et al. 2003 ; 2004), are closer to the Sun and in transition from 
the Kuiper belt towards the inner solar system. More than 1300 
TNOs have been detected so far, revealing a rich orbital structure 
and intriguing physical properties. The Trans-Neptunian pop- 
ulation is analogous to the debris disks observed around sev- 
eral other, 5-500 Myr old stars (Moro-Martin et al. 120081 Jewitt 
et al. 120091 1. This analogy is bolstered by similarities in sizes 
and observed masses (typically 30-300 AU and 0.01-0.1 M e for 
the "exo-disks"), with the important difference that the detected 
mass in extra-solar debris disks is in the form of ~ 10-1 000 jjm, 
short-lived, dust particles. The vast majority of the mass in these 
disks is invisible to us, probably in the form of kilometre (or 
more)-sized bodies, resembling trans-Neptunian objects. 

As part of the Herschel (Pilbratt et al. 2010) Science 
Demonstration Phase we observed 17 targets in different 
instrument configurations and observing modes. Here we 
present the analysis of early photometric measurements with 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



the Photodetector Array Camera and Spectrometer (PACS - 
Poglitsch et al. 1201 Oi l of five TNOs and two Centaurs. The sci- 
ence aspects from longer wavelengths Spectral and Photometric 
Imaging Receiver (SPIRE - Griffin et al. 1201 Oi l observations on 
(136472) Makemake and (90482) Orcus are included in Lim et 
al. ([2010) and the thermal lightcurve of (136108) Haumea is 
presented by Lellouch et al. ( 120101 ). The full Open Time Key 
Programme includes about 140 TNOs (Muller et al. 120091) . 



2. Observations and data reduction 

The PACS photometric measurements (70/100/160/im bands) 
were either taken in point-source mode with chopping-nodding 
on three dither positions (pre-launch recommended mode for 
point-sources), or in mini scan-map mode covering homoge- 
neously a field of roughly 1' in diameter (Poglitsch et al. 120101 . 
The mini scan-map mode turned out to be more sensitive and 
better suited for our project. 

The chop-nod data reduction was done in a standard way 
(Poglitsch et al. 1201 01) . The scan map processing deviated from 
the default way: Two scans were joined (for combined scan- 
maps) before executing first an unmasked high-pass filtering to 
identify sources and bright regions, which were then masked for 
a second high-pass filtering (without deglitching on the masked 
sources). The filter widths of the high pass were empirically cho- 
sen to be 62" and 82" in for the 70/100 //m maps and 160 jum 
maps, respectively, for S/N and flux conservation reasons. Then 
second order deglitching was applied in the source regions and 
the final maps were created. The calibration was done by apply- 
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Table 1. Observation summary: target name, observation identifier from the Herschel Science Archive (OBSID), observation mid- 
time (UT) in 2009, observation duration [s], mode: chopping/nodding/dithering or map parameters for 20"/s mini scan-map mode 
(scan leg length, separation, number of scans, map orientation in detector array-coordinates), repetition factor either for a chop-nod 
cycle or for full map, PACS photometer bands, colour-corrected flux values at PACS photometer reference wavelengths A c : 70 or 
100, as well as 160 /mi. 



Target 



Obs. ID Mid- time [UT] Dur. [s] Mode/Remarks 



Bands 



FD [mJy] 



(208996) 2003 AZ84 


1342187054 


11-16 19:20 


2526 


chop/nod/dither/ 1 6 


70/160 


27.0±2.7 


19.8+5.2 


(126154) 2001 YH140 


1342187062 


11-17 18:37 


5666 


chop/nod/dither/3 6 


70/160 


9.8+2.9 


<13 


(79360) 1997 CS29 


1342187073 


11-18 14:24 


5666 


chop/nod/dither/3 6 


70/160 


5.1±1.2 


14.5+2.9 


(82075) 2000 YW134 


1342187074 


11-18 16:00 


5666 


chop/nod/dither/3 6 


70/160 


<5 


<8 


(42355) Typhon 


1342187113 


11-20 00:05 


1584 


chop/nod/dither/ 1 


70/160 


17.0±3.4 


<13 




1342187114 


11-20 00:33 


1584 


chop/nod/dither/ 1 


100/160 


16.4+1.9 


<10 




7113&7114 


11-20 00:19 


3168 


combined 


160 




<9 


2006 SX368 


1342188416 


12-21 19:02 


2722 


3.5'/4.0"/10/ 63°/9 


70/160 


24.1 + 1.3 


<15 




1342188417 


12-21 19:49 


2722 


3.574.0710/11779 


70/160 


26.4+5.8 


<22 




8416&8417 


12-21 19:26 


5444 


combined scan-maps 


70/160 


22.2+2.9 


<17 


(145480) 2005 TBI 90 


1342188482 


12-23 20:08 


2722 


3.574.0710/ 6379 


100/160 


4.6+0.7 


<7 




1342188483 


12-23 20:55 


2722 


3.574.0710/11779 


100/160 


5.5+0.8 


3.6+1.5 




8482&8483 


12-23 20:24 


5444 


combined scan-maps 


100/160 


4.7+0.6 


<6 



ing flux overestimation corrections of 1.05, 1.09 and 1.29 at 70, 
100 and 160//m, as recommended in the PACS release note on 
point-source photometrjQ. 

In order to obtain monochromatic fluxes at the reference 
wavelengths 70, 100 and 160/mi we applied colour corrections 
of 0.98, 0.99 and 1.01 (Poglitsch et al. l20T0l l. with uncertain- 
ties of approximately + 1 -2% related to the full range of possible 
TNO and Centaur colour temperatures. The photometry of the 
targets was done by applying a standard technique, described 
in Lellouch et al. (120 1 Oi l. Table Q] summarises the selected SDP 
PACS observations with relevant instrument and satellite pa- 
rameters together with the obtained fluxes. The monochromatic 
fluxes are listed with lcr-errors, including also systematic errors 
like sky background gradients and variations. For non-detections 
we gave the 3cr errors as upper flux limits. 



3. Observational results and model input 
parameters 

In order to derive sizes and albedos from thermal-IR observa- 
tions a thermal model is required. The thermal emission of an 
atmosphereless spherical body at distance A is given by 



F(A) 



eR 2 /A 2 



II 



B[T(0,<p)]cos 2 <pcos(0-a)d0dip (1) 



where e is the emissivity, R the radius of the object, B the 
Planck function, ip the latitude, 9 the longitude measured from 
the sub-solar point, and a the solar phase angle. As is usual for 
small bodies we adopt e = 0.9 throughout this work. The use 
of this equation requires a model of the temperature distribu- 
tion T(0, <p) over the surface. We applied the following models: 
(i) The standard thermal model (STM, Lebofsky et al. 1986), 
which describes a non-rotation, low thermal inertia surface in 
instantaneous thermal equilibrium with insolation (beaming pa- 
rameter 77=0.756, optimised for mid-IR colour temperatures of 
large main-belt asteroids); (ii) The fast rotating model or isother- 
mal latitude model (FRM, ILM, Veeder et al. 1989 Lebofsky 
& Spencer 1989), as an extreme case of a high thermal inertia 
and/or fast rotating body; (iii) Intermediate models with respect 
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to STM and ILM, e.g., the near-Earth asteroid thermal model 
(NEATM, Harris [T9"9"8"j l, where the beaming parameter 77 is fitted 
to the multi-band data (T oc 77~'^ 4 ); (iv) A thermophysical model 
(TPM, Lagerros [T9961 [19971 [199%1 Miiller & Lagerros [19981 . 
where the temperature distribution is calculated for the given 
illumination and observing geometry, rotation axis and period 
have to be assumed if not available. The "free" parameter is the 
thermal inertia, a physical property of the surface material. The 
TPM assumes a "default" roughness (Miiller & Lagerros l2002l l. 
while the 77-driven models work with smooth surfaces. 

Table [2] summarises the observing geometries, the H- 
magnitudes with absolute uncertainties, lightcurve influences 
and rotation periods which were used as input for the modelling. 
Effective diameter D e fr, geometric albedo py and Hy magnitude 
are connected via D eff = 1347.4 x 10~ 02H 7 ^Jp^, (Bowell et al. 
1989). Knowing Hy (i.e., the reflected part of the Sun-light) and 
measuring fluxes at thermal wavelengths (i.e., the thermally re- 
emitted Sun-light) allows now to solve for D e ff and py, the ef- 
fective diameter of the target and its geometric albedo. An un- 
certainty of 0. 1 mag in Hy leads to errors of 4-5% in the derived 
diameter and 8-10% for the geometric albedo. 



4. Results and discussion 

The combined Herschel and Spitzer data show that neither of the 
"canonical" models (STM, ILM) fit the observed fluxes over the 
entire wavelength range (see Fig.[TJand Fig. [2]). The data require 
either an intermediate beaming value 77 (or different ones for dif- 
ferent parts of the SED) or a more sophisticated TPM with the 
thermal inertia F as key parameter for the temperature distribu- 
tion on the surface. The value 77 = 1 corresponds to a smooth 
surface with zero thermal inertia. For our observations, which 
are carried out at low solar phase angle, thermal inertia would be 
expected to raise 77 (because it reduces the day-side temperature), 
while roughness leads to higher-than-expected effective temper- 
atures, hence it lowers 77. Both model techniques were applied. 
For the TPM we used a ^-technique to find optimum solutions 
and uncertainties which are compliant with the observed fluxes 
and errors. The 77-NEATM uncertainties are either based on a 
bootstrap Monte-Carlo analysis (e.g., Mueller et al. 2010) for the 
three targets where MIPS and PACS observations are available 
or on a fixed 77 = 1.2 + 0.3 in all other cases. 
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Table 2. Herschel observing geometries for our seven targets, including properties derived from ground-based visible observations 
and the obtained radiometric solutions, r: Sun-target distance, A: Herschel-target distance, a: phase angle; Hy magnitudes, lightcurve 
A mag , and the rotation period P; references for the preceeding three columns, derived effective diameter Dj™ in [km] and geometric 
albedo p^ PM values from the TPM analysis, the corresponding possible thermal inertias T [J irT 
value (*: for a fixed value). 



s- ui K" 1 ], and the fitted NEATM 77 



Target 



r[AU] A[AU] a[°] Hv [mag] A„ 



P[h] 



Ref. 



r^TPM 



TPM 



2003 AZ 84 


45.376 


44.889 


1.11 


3.83±0.04 


0.14±0.03 


6.79 h 


1,2,3 


850-970 


0.05-0.09 


2-10 


1.31 


2001 YH 140 


36.576 


36.169 


1.44 


5.8±0.2 


0.13+0.05 


13.2h 


4 


300-390 


0.06-0.10 


0-10 


1.2* 


1997 CS 29 


43.509 


43.241 


1.27 


5.65±0.09 


<0.08, <0.22 


— 


5, 6, 7, 8 


250-420 


0.06-0.14 


0-25 


1.2* 


2000 YW I34 


44.125 


43.833 


1.25 


4.88±0.05 


<0.10 


— 


9, 10 


<500 


>0.08 


0-25 


— 


Typhon 


17.943 


18.376 


2.84 


7.68±0.04 


0.07+0.01 


9.67 h 


4, 11, 12 


134-154 


0.065-0.085 


1-10 


0.96 


2006 SX 168 


11.972 


12.271 


4.48 


9.5 


— 


— 


MPC 


70-80 


0.05-0.06 


0-40 


1.2* 


2005 TB 190 


46.377 


46.683 


1.17 


4.58±0.22 


0.12+0.02 


12.68 h 


13 


335-410 


0.15-0.24 


0-25 


1.2* 



Refere nces: (1) Fornasier et al. [20041 ( 2) D eMeo et al. 2009] (3) Pe rna et al. |2010; (4) T hirouin et al. l20T0l (5) Shepp ard & 

Jewitt 120021 (6) Romanishin & T egler [T9991 (7) Davies e t al. [20001 (8) Boehnhard t et a l. 120011 (9) Sheppard & Jewitt l2"003l 

(10) Doressoundiram et al. 2005, (11) Rabinowitz et al. 2007; (12) Tegler et al. 2003; (13) Thirouin et al. in prep.; MPC: 
http : //www . cf a . harvard . edu/iau/lists/Centaurs . html 
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Fig. 1. Observed Herschel and Spitzer flux values for 208996 
(2003 AZ84). This example illustrates that the simple "canoni- 
cal" STM and FRM/ILM fail to match the full SED range. Either 
a model with floating r\ or a TPM is necessary for the observed 
TNOs. 



(208996) 2003 AZ 84 is a Plutino in 3:2 mean-motion reso- 
nance (MMR) with Neptune with an eccentricity^ of 0.18 and 
inclination of 13.5°. In combination with fluxes from Table [TJ 
we used updated 24 and 70/im MIPS fluxes from Stansberry 
et al. (2008) (F 24 =0.28±0.02mJy; F 70 =24.6±3.1 mJy). The best 
77-based model solution resulted in an effective diameter of 
896±55km and a geometric albedo of py=0. 065+0. 008, the cor- 
responding beaming parameter is 7/=1.31±0.08. A more sophis- 
ticated TPM analysis on the basis of a spherical body with the 
given rotation period and a spin axis perpendicular to the solar 
direction, supported by a strong visible lightcurve (see Table 0, 
gave very similar D e ff and pv values at a thermal inertia of 
5^ J rrT 2 s~ 05 K" 1 . The very low thermal inertia is an indication 
that the surface might be covered by loose regolith with a low 
heat capacity in poor thermal contact. Although there are indica- 
tions of crystalline water ice on the surface (Guilbert et al. 2009), 
a solid compact layer of ice can be excluded, because this would 
require significantly higher thermal inertia. It was also possible, 



via the TPM, to investigate the influence of spin axis orienta- 
tion, sense of rotation and rotation period, but for objects with 
very low thermal inertia the influence of the rotational proper- 
ties are only on a level of less than 5% on the effective diameter 
and geometric albedo solutions. The diameter and albedo results 
agree with the Stansberry et al. ([2008) values within the given 
error bars. The results of all model fits are shown in Fig.Q] The 
lightcurve influence (see A mag in Table |2]l is significant and has 
been taken into account (0. 14 mag correspond to a ~7% diameter 
change and ~14% albedo change). 



(126154)2001 YH140 is a dynamically hot object with a semi- 
major axis close to the 5:3 MMR with Neptune, which might 
have excited the orbit from a previously dynamically colder or- 
bit. 2001 YH140 was not observed by Spitzer and our PACS ob- 
servation was done in chop-nod technique and with only a single 
epoch in the 160 pm band. Based on an assumed 77 = 1.2+0.3 the 
70 pm detection leads to a NEATM-solution of D efl -=349+81 km 
and /?y=0.08+0.05. The result from the TPM analysis is similar 
with a possible diameter range between 300 and 400 km and ge- 
ometric albedos between 0.05 and 0. 10, based on a range of ther- 
mal inertias from to 10 J irT 2 s~°- 5 K . These solutions corre- 
spond to predicted flux values of 6-1 1 mJy at 160/im, compati- 
ble with the upper flux limit of 13 mJy at 160/im. 



orbit parameters as FK5/J2000.0 helio. ecliptic osc. elements 



(79360) 1 997 CS29 is a dynamically cold object with at least 
one satellite (Stephens & Noll 2006) and with an orbit very 
close to the Neptunian 7:4 MMR. The NEATM analysis of the 
two PACS measurements was not conclusive, only an unreal- 
istically high 77-value above 5 could connect the two fluxes. 
But this solution would correspond to an object with more 
than 1000km diameter and an albedo below 1%. Combining 
the PACS 70pm observation with two Spitzer measurements 
(F 24 =0.057±0.005mJy, F 70 =3.0+01. 6 mJy) confirmed that the 
160 pm flux is very likely contaminated from an unknown back- 
ground source. The radiometric analysis via TPM techniques 
(excluding now the 160/im flux) leads to a possible diameter 
range between 250 and 420 km and an albedo of 6-14% for a 
large range of thermal inertias. The TPM analysis of the Spitzer 
observations alone are indicative of a very low thermal inertia 
at around 1 JirT 2 s~ a5 Kr', for a D e ff = 270 km and py=0.14. 
Here, the fixed r\ approach gave 402±69km and 0.06±0.02, 
favouring the higher inertias above 10 J m~ 2 s~ 05 K" 1 . 
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(82075) 2000 YW134 is a binary which moves on an extremely 
eccentric orbit and is dynamically a detached object rather than 
a member of the scattered disk. In a 1.5 h measurement we only 
obtained upper flux limits. Nevertheless, the observation (mainly 
the JOjjm flux) combined with the Hy magnitude in Table[2]con- 
strain the possible diameter and albedo solutions: The target has 
to be smaller than about 500 km and the geometric albedo higher 
than 8%. Lower albedos or larger diameters would not be com- 
patible with the non-detection at 70 //m. 
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Fig. 2. The observed Herschel and Spitzer flux values for 
Typhon. The line-styles are the same as in Fig.Q] 



(42355) Typhon is a binary Centaur on a highly eccentric 
orbit influenced by Uranus (Alvarez-Candal et al. 20101. We 
were using all Spitzer and PACS detections to derive radio- 
metric properties. Via NEATM we obtained D e ff=138±9km, 
Pv=0.080+0.01 for a 77-value of 0.96+0.08. At a thermal inertia 
of 5^ J irT 2 s~° 5 K the measured fluxes and the corresponding 
TPM predictions agree very well (Fig. [2]). The derived effective 
diameter is 144±10km, the geometric albedo py=0.075±0.010. 
The 3<x upper flux limits at 160^m indicate that the diame- 
ter must be around 140 km or smaller, the albedo at 0.08 or 
larger, both values agree with the derived radiometric properties. 
Larger thermal inertias can be excluded due to the non-detection 
at 160/im, very small thermal inertias would cause difficulties 
to match both Spitzer fluxes. Combined with the binary system 
mass derived by Grundy et al. (120081) . we derived a bulk density 
of 0.66+j] Qg g cirT 3 , which is slightly higher than the Grundy et 
al. (2008) value which was based on a single-band Spitzer detec- 
tion (Stansberry et al. 120081 ). 

2006 SX368 is a Centaur on a very eccentric orbit, near the 5:4 
MMR with Uranus. The PACS 70//m detection is compatible 
with diameters in the range of 70-80km and an albedo of 0.05- 
0.06, allowing for a wide range of thermal inertia between 0- 
40 JirT 2 s~ - 5 K . The upper 160yum flux limit constrains the 
diameter range to values below 105 km and an albedo larger than 
0.03, both in agreement with the measured 70|/m-flux. The fixed 
77-approach gave D e ff=79±9km and pv=0.05±0.01. 

(145480) 2005 TB190 has a highly eccentric orbit and belongs 
to the detached objects. Its aphelion lies beyond 106 AU, where 
it may even pass through the termination shock and into the he- 
liosheath on each orbit. Both PACS bands are beyond the emis- 
sion peak for (145480) 2005 TBi 90 and constrain the TPM out- 
put to an effective diameter in the range 335-410km and an 



albedo in the range 0.15-0.24. The ~0.3mag uncertainty in Hy 
combined with lightcurve variations does not influence the diam- 
eter solution significantly, but extends the possible albedo range 
to 0.12-0.30. The fixed 77-approach gave D e fF=375±45km and 
py=0.19±0.05. It is the highest albedo object in our sample. 

5. Conclusions 

Our small and dynamically very inhomogeneous sample con- 
firms the consistency between different model techniques and 
nicely agrees with Spitzer results (Stansberry et al. 2008]) on 
three overlap targets. Based on the seven targets we also showed 
the model capabilities for multiple, dual, single or even non- 
detections. Models with either a beaming parameter of ~ 1.2 or 
thermal inertias below 25 JirT 2 s~ - 5 K _1 explain the measured 
SEDs, confirming low heat conductivities at temperatures far 
away from the Sun. The target sizes range from diameters be- 
low 100 km to almost 1000 km. The derived geometric albedos 
are below 10%. Only (145480) 2005 TBi 90 has an albedo above 
15%, which is possibly related to its very unusal orbit located 
entirely outside the major planets. The derived low albedos are 
close to the mean value of 8% given by Stansberry et al. (|2008). 
Albedo trends with object size, dynamic or taxonomic types are 
not yet visible in the small sample. 
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